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Abstract - The instanton-noninstanton (I-NI) transition in the tunneling process, which has 
been numerically observed in classically nonintegrable quantum maps, can be described by a 
perturbation theory based on an integrable Hamiltonian renormalized so as to incorporate the 
integrable part of the map. The renormalized perturbation theory is successfully applied to the two 
quantum maps, the Henon and standard maps. In spite of different nature of tunneling in the two 
systems, the I-NI transition exhibits very common characteristics. In particular, the manifestation 
of I-NI transition is obviously explained by a remarkable quenching of the renormalized transition 
matrix element. The enhancement of tunneling probability after the transition can be understood 
as a sudden change of the tunneling mechanism from the instanton to quite a different mechanism 
supported by classical flows just outside of the stable-unstable manifolds of the saddle on the top 
of the potential barrier. 


Poincare proved that almost all the classical Hamilto¬ 
nian systems are nonintegrable. The tunneling process 
in classically integrable systems is described almost com¬ 
pletely by the instanton theory and is understood in terms 
of classical trajectories [T|. However, the theory of tunnel¬ 
ing in classically nonintegrable systems remains far from 
complete, and is still one of the unsolved fundamental 
problems in theoretical physics. It has been more than 
two decades since the study of nonintegrable tunneling 
starts EH z] ( see recent progress in 00). In the energy 
domain approach, the application of the Herring’s formula 
to nonintegrable systems mam], hybrid approaches com¬ 
bining pure quantum and semiclassical theory, which was 
applied to explaining chaos-assisted tunneling mmm 
and resonance-assisted tunneling ehui, and so on have 
been proposed. On the other hand, in the time domain 
approach, the complex domain semiclassical theory has 
revealed a crucial role of the complexified stable-unstable 
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manifold mechanism, which implies the importance of Ju¬ 
lia sets U5HHI- However, the relation and correspondence 
between theories and numerical investigations are not yet 
clear enough and there is no unified view. 

Most of works have been devoted to the so-called quan¬ 
tum maps which model the essence of nonintegrable fea¬ 
ture of quantum dynamics in a simple minded way. A 
very important result which has been reported in several 
works is that even a very weak nonintegrable perturba¬ 
tion leads to a remarkable enhancement of the tunneling 
rate 

mmmuBHS]- This transition is a common feature 
observed for quantum maps but its basic origin is still 
controversial and is not understood completely. 

Wave functions of non-integrable systems are localized 
on classical tori in the nearly-integrable regime according 
to the well-known Einstein-Brillouin-Keller (EBK) quanti¬ 
zation rule. A remarkable feature of the invariant torus in 
nonintegrable system is that it breaks at a border called 
the natural boundary on which singularities are densely 
accumulated [22] , 

It was conjectured that the natural boundary interrupts 
the instanton, thereby may influence the tunneling rate 
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0 . Indeed, clear evidence has been demonstrated recently 
manifesting that the remarkable enhancement of the tun¬ 
neling rate is closely connected with the interruption of 
instanton due to the natural boundary m- We call the 
first transition from the instanton tunneling to some un¬ 
known type of tunneling as the instanton-noninstanton (I- 
NI) transition, which occurs in a notable manner in quan¬ 
tum maps. In the present paper, by introducing a pertur¬ 
bation theory based upon maximal renormalization of the 
integrable part, we show clearly that a different tunnel¬ 
ing process supported by classical flow just outside of the 
stable-unstable manifolds of the saddle at the top of the 
potential barrier induces the transition to the NI regime. 

In what follows, we take the symmetric form of the 
quantum map U = e ~ tp2 / 4h e ~ leV (Q)/ h e - tp2 / 4h ( 0 r U = 
e -ieV(Q)/ 2 h e -iP 2 / 2 h e -ieV(Q)/ 2 hy Redefining the effective 
Planck constant as n = and the new set of conju¬ 

gate operators p = -ind/dq and q = Q 1 then U becomes a 
convenient form in our approach 

U - e ~iyAp 2 /4:K. e -iy/eV(q)/K, e -iy/ep 2 /4K' 

Figs. Oja) and (b) are typical examples of tunneling 
characteristics computed for two typical quantum maps, 
namely the Henon and standard maps. The figures show 
the representative tunneling probability as a function of 
the quantum number, and they exhibit a typical feature 
of the I-NI transition. Here the quantum maps have the 
potential functions V(Q) = 2 Q 2 + Q 3 /3 (the Henon map) 
and cos Q (the standard map), respectively. Although 
these examples are nonintegrable, we here examine weakly 
nonintegrable regimes with small e as displayed in the in¬ 
sets of fig. [Q Therefore, the classical invariant tori still 
remain according to the celebrated Kolmogorov-Anold- 
Moser (KAM) theory, and good quantum numbers can be 
assigned to each quantum eigenstate following the EBK 
quantization rule. 

In the limit of vanishing nonintegrability, both systems 
are approximated by one-dimensional systems with poten¬ 
tial barriers which classically bound the eigenstates but 
quantum mechanically allow to tunnel by the instanton 
mechanism as is usual in one-dimensional systems, which 
will be discussed later. However, as the quantum number 
n decreases from the highest exited state with the energy 
just below the potential barrier height, a transition occurs 
at a characteristic quantum number, denoted hereafter by 
n c, and the tunneling probability suddenly deviates from 
the instanton probability forming a plateau. The com¬ 
mon feature seems to be quite paradoxical in the sense 
that the low lying states distant from the saddle point lo¬ 
cated on the barrier top, which is the very origin of the 
classical nonintegrability exhibit non-instanton tunneling, 
while the highly excited states close to the saddle point 
obey the instanton tunneling. Moreover, in both exam¬ 
ples the number of the eigenstates in the instanton regime 
n m ax - n c is insensitive to h , where n max is the maximum 
quantum number of the classically bounded states inside 



Fig. 1: The tunneling probability (in log 10 ) for the eigenstates 
of eq. (1) vs the quantum number n (less than n max ) for (a) 
the Henon and (b) the standard map with various values of h 
in the unit of ho and the nonlinear parameter t. The tunneling 
probability is computed in the asymptotic region (Q « -4) 
for the Henon map and at the the potential top (Q = 0) for 
the standard map (only the even parity eigenfunctions in the 
^-representation are considered). Here ho = 0.63>/e for the 
Henon map and ho = || s/e for the standard map. Insets shows 
the phase space portrait of classical maps with e = 0.1 and 
e = 0.36 for the Henon and standard maps, respectively. 


the potential well. 

To understand the I-NI transition the most natural 
way is to apply a perturbation theory based on an in¬ 
tegrable limit which has the instanton as the tunneling 
mechanism. The crudest integrable approximation of sym¬ 
metrized U in the small limit of the nonlinear parameter 
e is Hi = where Hi = p 2 /2 + V(q), but the dif¬ 

ference \U - Ui\ ~ e 2 /n is too large to afford any signif¬ 
icant result for the exponentially small tunneling effect. 
We therefore develop a systematic expansion with respect 
to the smallness parameter e which renormalizes the in¬ 
tegrable part of U into a single effective Hamiltonian as 
much as possible. A possible candidate to achieve this 
is the Baker-Hausdorff-Campbell (BHC) expansion, which 
approximates the product of exponential operators m sys¬ 
tematically in terms of a single exponential operator ex¬ 
pressed by an effective Hamiltonian, which coincides with 
Hi in the lowest-order approximation. 

Below, we apply a renormalized perturbation theory to 
the Henon and standard maps. They are both simplest 
classes of maps exhibiting quite different nature of tunnel¬ 
ing, which is understood by the simplest one-dimensional 
Hamiltonian Hi = p 2 /2 + V(q): the Henon map has the 
cubic potential with the bottom at q = 0 and tunneling 
is an irreversible transport of probability toward q = -oo 
going over the barrier at q - -4. For the standard map the 
47 t periodic boundary condition is imposed in the q direc¬ 
tion. Then its cosine potential gives two symmetric valleys 
with bottoms at q = ±n and separated by two symmetric 
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barriers at q = 0 and q = 2n. 

We can show that the BHC expansion of U leads to an 
effective polynomial Hamiltonian 

fltf- E (2) 

ke,iek,jik 

with H e = Yj K k C(i, k, i,j)p l q\ 
k,i,j =0 


where C(n,k,i,j) is the set of coefficients of 0(1), and 
the terms of k - 0 provide the classical Hamiltonian. Note 
that only the odd powers of e 1 / 2 appear in the sum because 
of the symmetrized form of U. The unitary operator thus 
induced is given as Um = e . We do not intro¬ 

duce any artificial absorbers and/or absorbing boundary 
conditions, which may crucially influence the original dy¬ 
namics mm. 

Once the renormalized Hamiltonian is obtained, we 
straightforwardly develop a perturbation theory taking 
the difference A Um = U - Um = 0{e ( ‘ M+l ' > sfe/ k) as the 
perturbation. Then the lowest-order perturbative eigen¬ 
function is given as 




k M ViA< M) ), 

„ (u^AUmIu^) 

i e -V^i M) _ e -V^r ) 


,( M )\ 


(3) 


where \u^ I ' > ) and E„ M ^ are respectively the eigenfunction 
and the energy eigenvalue of the ID integrable Hamilto- 
man . 

If e is not so large, there appear two sorts of fixed points 
in classical phase space: stable and unstable fixed point 
(0,0) and (-4,0) for the Henon map, and (±7r,0) and 
(0,0) for the standard map respectively. As shown in 
the insets of fig. [T] KAM tori predominate phase space 
in both cases and KAM regions are encircled by the sta¬ 
ble manifold W s and the unstable manifold W u of the 
saddle fixed point. KAM tori support quantum eigen¬ 
states, each of which satisfies the EBK quantization con¬ 
dition f pdq/2r: = (n + 1/2 )n (n = 0, l,-",n max )- They 
all have finite tunneling life-times (Henon map) or tun¬ 
neling splittings (standard map) due to the tunneling via 
the instanton trajectory, which is very well approximated 
by the lowest order instanton solution p = yj2(V(q) - E n ) 
or H\{\p,q) = E n ( E n is the quantized energy). A no¬ 
table feature of the renormalized Hamiltonian (p ; q) 

is that tunneling tails of eigenfunctions (q\u^^) do not 
change very significantly even if one increases the order of 
renormalization M. However, according to eq.dU), the ap¬ 
plication of perturbation changes drastically the tunneling 
tail of eigenstates with increase in M. 

In fig. [2jb) we show a typical result of renormalized per¬ 
turbation theory. For lower-order approximation, U - Um 
is too large to control correctly the exponentially small 
tunneling component and so the perturbative solution 


yields quite absurd results. But as the order M of renor¬ 
malization increases the perturbative solution converges 
to the exact eigenfunction |'3/ n ) obtained by numerical di- 
agonalization and reproduces even the complicated oscil¬ 
lations at the tunneling tail as demonstrated in figs. [2] 
(al) and (a2). The higher-order renormalized perturba¬ 
tion calculation also succeeded in reproducing the exact 
tunneling probability in a rather wide regime including 
both I (instanton) and I-NI transition regions. (See figs. 
[2] (bl) and (b2)). 



Fig. 2: The convergence of tunneling tail to the exact numerical 
one (black) with increase in M in case of (al) n = 13-th excited 
state of the Henon map with e = 0.1 and h = ho/2, and (a2) the 
ground state of standard map with e = 0.36 and h = ho/2. The 
exact and higher order renormalized perturbative results are 
not distinguishable with each other. Right-hand panels show 
the probability amplitude of tunneling tail of exact numerical 
diagonalization (black lines) and of perturbation calculation as 
a function of quantum number (less than n max ) in case of (bl) 
the Henon map with the order M = 21, and (b2) the standard 
map with the order M = 7. In (b2) the quantum number is 
shifted by -5 for ho, and so n = -5 means n = 0 for ho- 


The success of the perturbation theory means that the 
origin of the I-NI transition may be resolved into inte¬ 
grable bases: according to eq. © we define the contri¬ 
bution spectrum representing the amount of contribution 
from the j-th eigenstate of the integrable model to the 
n-th perturbative eigenfunction: 

Con^M(l(9l^ M) )l))l( u 5 M) l^ M) >l" 

|((l(9|w5 M) )l))l<^l A ^« M) >l ( if 3 * n ) (4) 

t<(l(9|wi M) }|)) (if j = n) 


where ((...)) means to take average over a range of q in the 
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asymptotic region (Henon map) or around q = 0 (standard 
map). 

We discuss closely how the renormalized perturbation 
theory describes the I-NI transition by taking the contri¬ 
bution spectrum of the Henon map as an example. Figure 
[3](a) shows how the contribution spectrum Coi/^2 (0 < 
j ^ n m ax) varies with the order of renormalization M at 
two representative quantum numbers n before and after 
the I-NI transition. A quite interesting fact is that for 
M = 1 the largest contribution comes from a broad peak 
centered at j = n*, where n* is the number of a quantum 
state just above the threshold of dissociation n = n max . 
However, with increase in the renormalization order M, 
the transition matrix elements \(u ^ M ^] AC/ )| are re¬ 
markably quenched, and the peak at j = n* is reduced in 
the logarithmic scale and finally overwhelmed by the very 
sharp peak at j = n, which means that the contribution 
from instanton of the integrable basis ( q\ui^) is predomi¬ 
nant in the tunneling process. On the other hand, at n less 
than a critical quantum number n c the quenching of the 
transition matrix element can no longer reduce the peak 
at j = n* less than the instanton peak, and the contribu¬ 
tion to tunneling is attributed to the broad peak around 
j = n*. The competition between the two peaks explains 
the characteristics of I-NI transition. It should be empha¬ 
sized that without remarkable quenching of the transition 
matrix element by renormalization the instanton phase is 
absent and the I-NI transition cannot be observed. 



Fig. 3: Typical contribution spectra Conj^ of (al) before and 
(a2) after the I-NI transition, (al) and (a2) respectively cor¬ 
respond to n = 15 and n = 7 of the Henon map with e = 0.1 
and h = ho/2 (see fig. [2](bl)). The renormalization order grows 
as M = 1,5,9,13,17,21. (b) The major transition matrix ele¬ 
ment \{u^^\AUm\u^' > )\ vs n for the same Henon map, where 
n* = 21 = n max + 2. M is increased from 1 to 21. The black 
doted curve denotes the instanton peak amplitude ((|(9|wi M ^)|)) 
of the n-th state. The vertical axis is drawn in the log 10 scale. 


Since the above threshold state \u^P) is connected with 
q = oo by real classical paths, |(g|u£^)| ~ 0(1) and so 
the peak strength Con^2 n at n* is approximated by the 
transition matrix element |(ul M ^|AI7)l- Figure [3] 


(b) plots \{ui M) \AU M \u^)\ at just above the threshold 
for various values of M. The instanton peak strength 
((|(g|u^ M ^ )|)) is also shown as a function of n, which does 
not significantly depend on M as mentioned above. 

With increase in M, the curve \(u^ I ' , \AUm\uI^ I ' > )\ vs n 
is largely deformed to show a very characteristic struc¬ 
ture: it decreases steeply as n decreases below n max , 
but it reaches a definite plateau and then it decreases 
again steeply. It is just on the plateau that the curve 
Ku^IAUmIu^I vs n intersects with the instanton peak 
curve, which means that tunneling using the transition to 
the state n = n* overwhelms the instanton tunneling and 
so the intersection determines n c . As M increases, the 
height of plateau decreases rapidly and reaches finally to 
a limit, which causes the instanton region n c < n < n max to 
grow from a null region to a finite region with the width 
proportional to \/e. Below n c , the transition matrix ele¬ 
ment controls tunneling and so the tunneling probability 
follows the plateau structure, which explains the charac¬ 
teristic plateau of the tunneling amplitude seen in fig[2] 

The origin of the transition to a plateau-like characteris¬ 
tics of the tunneling rate from the instanton tunneling rate 
reported in preceding works can therefore be attributed to 
the formation of the plateau of the transition matrix ele¬ 
ment and a drastic decrease of plateau height in higher- 
order renormalization. The I-NI transition in the standard 
map in hg[ljb) follows the same scenario [ 23 . 

We note that if one shifts each curve in Fig. [T| hor¬ 
izontally such that the maximal quantum number n max 
for each curve coincides with each other, they all show 
very similar characteristics, and are insensitive to the ef¬ 
fective Planck constant k. Such a feature can hardly be 
explained by classical objects e.g., nonlinear classical res¬ 
onances, which will be discussed in detail in our forthcom¬ 
ing papers [23] . 

Finally, we elucidate the classical dynamical significance 
of the plateau region which is characteristic after the TNI 
transition. The above analyses tell us that the states com¬ 
posed of the self-component j = n and the ones forming 
the broad peak around j = n* play as the principal compo¬ 
nent. We define here the principal component contributing 
to tunneling (PCT) and observe the phase space Husimi- 
plot to investigate the classical interpretation for the I-NI 
transition. The PCT is defined as the projection of 
onto the principally contributing subspace constructed as 
follows. Let } be eigenfunctions rearranged in de¬ 
scending order of the magnitude Corn^n in the contribu¬ 
tion spectrum, and consider the Af-dimensional subspace 
spanned by |i^ M ^) (1 < j < K). Let be the minimum 
number of K which makes the relative distance between 
the vector and its projection onto the above intro¬ 

duced subspace less than a small enough threshold value 
Tth (« 1), namely, the minimal K = K min such that 

((- iu„W', -» < 
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Fig. 4: Husimi plots of the PTC for the eigenfunctions of the 
Henon map with e = 0.10 and h = ho/2 (see fig. [Hbl)). (a) 
before (n = 15), and (b) after (n = 13) the I-NI transition, and 
(c) is the eigenstate at the edge of plateau (n = 7). (a) and 
(c) correspond to fig. EJal) and (a2), respectively. Contours 
are plotted in log scale, and the shaded regions indicates the 
highest levels. The red curve is the W s -W u complex of the 
saddle S (red circle), and the green dot E denotes the stable 
fixed point. The Husimi levels of (a) coincide entirely with 
the eigenfunction of \u[^ r ' > ) of H^\ which has its major body 
inside of the W s -W u complex, while in (b) and (c) the major 
body winds around W a , passes nearly to S and is blown to -oo 
following the W u . The maximum level of contour is 10 0 ' 5 in 
(a), suddenly drops to 1(T 4 ' J in (b) after the I-NI transition 


Then the space spanned by \u^ M ^) (1 < j < K mm ) consti¬ 
tutes the principally contributing subspace to tunneling. 
We take rth = 0.2, for example. We have to emphasize that 
M must be taken sufficiently large (M > 13 in practice) in 
order that the PCT is significant. 

The Husimi plots of the PCT before and after the TNI 
transition are depicted in fig. [4] The shadowed region 
indicates the region with the highest probability level. In 
the instanton regime the PCT coincides almost with the 
eigenfunction of the unperturbed integrable Hamil¬ 

tonian, and the Husimi plot of PCT in fig. |4ja) indeed 
traces the classical invariant circle shown as the shadowed 
region, and the tunneling component is represented by 
the monotonously decaying contours encircling the shaded 
quantized invariant circle. In particular, when observed in 
the g-coordinate, the tunneling tail is the line p = 0 pass¬ 
ing across the saddle S , which is nothing more than the 
instanton. 

On the other hand, fig. HKb) indicates that a drastic 
change occurs in PCT when the transition to NI happens: 
the main component, characterized by almost the same 
probability levels, runs along the stable manifold W s , at¬ 
tracted and repelled by the saddle S, finally runs away to¬ 
ward q = -oo along the unstable manifold W u , which man¬ 
ifests that the PCT contributing to the tunneling tail rep¬ 
resents classical flows just outside of the W s -W u complex. 
The PCTs for all the eigenstates forming the plateau take 
almost similar patterns but they approach more closely 
to the W s -W u complex, as the quantum number shifts to 
the edge of the plateau. Figure [41(c) shows that just at 
the eigenstate on the edge of the plateau its PCT coin¬ 


cides with the W s -W u complex. As is seen in fig. [TJ the 
tunneling amplitude decreases rapidly when n is less than 
the plateau edge quantum number. All the above features 
are common in the Henon and standard maps, which are 
more completely discussed in [23] . 

In the present paper we have developed a perturbation 
theory for nearly integrable quantum maps. This is based 
upon an integrable Hamiltonian which is constructed by 
maximally renormalizing the integrable part of the map. 
This was successfully applied to investigate the TNI tran¬ 
sition commonly observed in nearly integrable quantum 
maps such as the Henon and standard maps. A remarkable 
quenching of the highly renormalized transition matrix el¬ 
ements explains the origin of the abrupt change of tunnel¬ 
ing characteristics at the TNI transition. The PCT anal¬ 
ysis reveals that the tunneling mechanism in the plateau 
region after the TNI transition is due to the classical flow 
outside of the W s -W u complex is responsible for the tun¬ 
neling process after the transition and the flow coincides 
with the W s -W u complex at the edge of the plateau. This 
suggests the crucial role of classical dynamics related to 
the W s -W u complex, and elucidating the relation to the 
complex stable-unstable manifolds mechanism based on 
the complex-domain semiclassics is strongly desired. 

* * * 
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